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SUMMARY 

A composite grid was generated for flows in turbines. It consisted of the 
C-grid (or 0-grid) in the immediate vicinity of the blade and the H-grid in the 
middle of the blade passage between the C-grids and in the upstream region. 

This new composite grid provides better smoothness, resolution, and orthogo- 
nality than any single grid for a typical turbine blade with a large camber and 
rounded leading and trailing edges. The C-H (or 0-H) composite grid has an 
unusual grid point that is connected to more than four neighboring nodes in two 
dimensions (more than six neighboring nodes in three dimensions). A finite- 
volume lower-upper (LU) implicit scheme to be used on this grid poses no 
problem and requires no special treatment because each interior cell of this 
composite grid has only four neighboring cells in two dimensions (six cells in 
three dimensions). The LU implicit scheme was demonstrated to be efficient and 
robust for external flows in a broad flow regime and can be easily applied to 
internal flows and extended from two to three dimensions. 


INTRODUCTION 
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lation an accurate solution for the complex fluid flow phenomena in a turbine. 

A smooth and nearly orthogonal H-grid can be generated if the turbine blades 
are thin and have only a slight camber and sharp leading and trailing edges. 

But such geometry is rare for a turbine because of the high work factor desired 
from the turbine and the mechanical factors to be considered in designing the 
blades, which usually operate at high temperature and speed (ref. 1). Turbine 
blades are often designed to have substantial thickness and camber and rounded 
leading and trailing edges. Geometries are more complicated in the radial- 
inflow turbine than in the axial turbine. A typical radial-inflow turbine 
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rotor, for Instance, turns the flow In 90° deflection In the meridional plane 
while simultaneously turning the flow In the blade-to-blade plane. Two repre- 
sentative turbine geometries are used to examine the grid generation for tur- 
bines: a radial-inflow turbine rotor (e.g., ref. 2) and an annular turbine 

stator (e.g., ref. 3). 

For a typical turbine blade no single grid offers satisfactory grid prop- 
erties In the entire rotor or stator passage. With the H-grld It Is difficult 
to obtain a good boundary-layer resolution around the leading edge. If very 
fine meshes are used to resolve this, the number of grid points In the upstream 
region becomes excessive. Thus mesh points would be wasted. With the standard 
C-grld (or 0-grld), for a typical turbine geometry, the grid becomes very 
skewed because the blades are highly cambered and the blade-to-blade grid lines 
curve sharply In the midpassage region (refs. 4 and 5). 

A two-dimensional 0/H patched grid used In a turbine cascade computation 
(ref. 6) shows the slope discontinuity at the 0-H grid Interface and has 
extremely large-aspect-ratio meshes near the leading edge. The numerical 
scheme used In that work Is a cell-centered scheme based on the Beam and Warm- 
ing approximate factorization (ref. 7). Although the flexibility of the scheme 
Is Indicated In two dimensions, Its likely limitations In three dimensions sug- 
gest an alternative approach. 

In this study two- and three-dimensional composite grids were generated 
In an attempt to Improve the grid quality In terms of smoothness, resolution, 
and orthogonality. A numerical scheme that will run on this composite grid Is 
discussed. 


STANDARD GRIDS FOR TURBINE GEOMETRIES 

A three-dimensional Invlscld analysis using a H-grld was performed for the 
radial-inflow turbine rotor (ref. 8). The rotor used In the analysis Is shown 
In figure 1. The rotor had 12 full blades thick enough to allow Internal cool- 
ing passages and tralllng-edge coolant ejection. The rotor had unswept radial 
blades at the Inducer Inlet and large blockages resulting from the thick trail- 
ing edges. In figure 2 the rotor surface Is defined by using some of the grid 
lines on the rotor surface. Three-dimensional meshes In the radial turbine 
passage are shown In figure 3. As an Illustration, only a selected number of 
grid lines on a selected number of grid surfaces are shown In the figure. A 
detailed H-grld for a blade-to-blade surface Is shown In figure 4(a). The mesh 
on a meridional plane Is shown In figure 5(a). Nonuniform grid spacing was 
used In all three directions. The quasl-orthogonal mesh lines were clustered 
near the leading and trailing edges, as Illustrated In figures 4(a) and 5(a). 

In addition, this feature was also used to position the Initial upstream grid 
at any desired location with a selected number of orthogonal grid lines and to 
position the final downstream grid where the experimental data were measured so 
that outputs from the Invlscld calculation could be directly compared with 
experimental results. In the H-grld, cusps were required at the leading and 
trailing edges to get convergence of flow solution. Figure 4(a) shows the grid 
skewness for the leading edge; figure 4(b) shows the rapid velocity change near 
the leading edge. In spite of the difficult blade geometry (l.e., thick edges, 
large trailing blockage, turning of the passage In the meridional and blade-to- 
blade directions simultaneously, etc.), the three-dimensional Invlscld analysis 
predicted the flow field reasonably well In a qualitative sense. But a large 
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entropy change occurred around the leading-edge region. As an example of 
results obtained from the three-dimensional Invlscld analysis, figure 5(b) 
shows secondary velocity vectors on cross-channel surfaces. 

A stator vane ring Is shown In figure 6(a). The full annular ring con- 
sisted of 36 vanes. A sector of four vanes used to make laser anemometer meas- 
urements Is shown. This vane ring was used for an annular cascade experiment 
reported In reference 3. The stator vane geometry Is shown In figure 6(b). 

The vanes were untwisted, had a constant profile from hub to tip, and had a 
height of 38.10 mm and an axial chord of 38.23 mm. The vane aspect ratio and 
the solidity at the mean radius (based on axial chord) were 1.0 and 0.93, 
respectively. The stator hub-tip radius ratio was 0.85 and the tip diameter 
was 508 mm. The stacking axis of the vane was located at the center of the 
tralllng-edge circle. To alleviate numerical errors associated with the H-grld 
skewness for the blunt leading edge, a standard C-grld was generated to examine 
any possible advantage over the H-grld for this turbine vane geometry 
(fig. 7(a)). The grid shows large skewing on the suction side of the blade 
passage (fig. 7(b)). Results of numerical simulations using this grid are 
likely to be affected by the grid skewing, as Indicated In references 5 and 6. 


COMPOSITE GRID 

Two- and three-dimensional composite grids were generated for the turbine 
vanes of reference 3 (figs. 8(a) and (b), respectively). The blades were 
untwisted and stacked at the trailing edge (fig. 8(c)). Only a selected number 
of grid lines are shown for Illustration. This composite grid consists of the 
C-grld (or 0-grld) In the Immediate vicinity of the blades and the H-grld In 
the upstream region and In the middle of the blade passage between the C-grlds. 
The C-grld (or 0-gr 1 d ) portion can be generated by using either the elliptic 
method (ref. 9) or the algebraic method (ref. 10). At the C-grld and H-grld 
Interface (or 0-H Interface) the slope continuity was preserved so that no 
special numerical approximations are needed for the derivatives at the Inter- 
face. The C-grld (or 0-grld) Is orthogonal to the blade surface and provides a 
good boundary-layer resolution near the leading edge. This composite grid has 
better smoothness, resolution, and orthogonality than any type of single grid 
for a typical turbine blade with large camber and rounded leading and trailing 
edges. 

For the two-dimensional grid In figure 8(a) the C-grld In the vicinity of 
the blade was generated by the elliptic grid generation code (ref. 11). Only a 
portion of the C-grld near the blade was retained for the composite grid. In 
the C-grld shown In figure 9, for Instance, the portion between n = 0 and 
n = n c < n max was retained, where n = 0 and n = n max are Inner and outer 
boundaries, respectively, and choice of nc Is arbitrary and depends on the 
extent of the shear flow region. Then the H-grld was smoothly patched with the 
C-grlds on the pressure and suction sides. At the C-H Interface the slope con- 
tinuity was preserved. The cubic spline and a stretching function were used to 
generate the H-grld In the midpassage and upstream regions. The three- 
dimensional grid In figure 8(b) was constructed algebraically from the two- 
dimensional composite grids of the hub and shroud surfaces. 
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NUMERICAL SCHEME 


Semidiscrete Finite-Volume Scheme 

The two-dimensional C-H composite grid has an unusual grid point that Is 
connected to more than four neighboring nodes (fig. 10(a)). At this point the 
usual differencing techniques cannot be applied. If a standard finite- 
difference scheme were used to solve a flow problem on this composite grid, It 
would be difficult to treat this special point. A finite-volume scheme 
(ref. 12) to be used with this composite grid presents no problem and requires 
no special treatment because each Interior cell of this composite grid has only 
four neighboring cells In two dimensions and six neighboring cells In three 
dimensions (fig. 10(b)). The finite-volume scheme Is described briefly here. 

The Euler equations In Integral form can be written as 

fff" in* Iff ■ as - o (i) 

n an 

for a fixed region £2 with boundary an. Here w represents the 
conserved quantity, F Is the corresponding flux term, and t Is time. 

A convenient way to ensure a steady-state solution Independent of the time 
step Is to separate the space and time discretization procedures. In the semi- 
discrete finite-volume scheme one begins by applying a semidiscretization In 
which only the spatial derivatives are approximated. To derive a semidiscrete 
model that can be used to treat complex geometric domains, the computational 
domain Is divided Into quadrilateral cells. Assuming that the dependent vari- 
ables are known at the center of each cell, a system of ordinary differential 
equations Is obtained by applying equation (1) separately to each cell. These 
have the form 


dt <S 1J w 1.i ) 


‘ ' 0 


where Sij Is the cell area and Qij Is the net flux out of the cell. 
This can be evaluated as 


( 2 ) 
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Z <ay k f k - 4* k g k ) (3) 

where f^ and g^ denote values of the flux vectors f and fl on the kth 
edge, Axk and Ayk are the Increments of x and y along the edge with 
appropriate signs, and the sum Is over the four sides of the cell. The flux 
vectors are evaluated by averaging the values In the cells on either side of 
the edge: 


f l -2 < f 1*l,3 * f 1,J> 

for example. The scheme constructed In this manner reduces to a central 
difference scheme on a Cartesian grid and Is second-order accurate In space 
provided that the mesh Is smooth enough. It also has the property that uniform 
flow Is an exact solution of the difference equations. 
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LU Implicit Scheme 


The composite grid Is generated to be used for both Invlscld and viscous 
flow calculations. For viscous calculations the C- or 0-meshes must be very 
fine to resolve the boundary layer. And It Is likely that the time step 
Imposed by an explicit stability bound will be much less than that Imposed by 
the accuracy bound of an Implicit scheme. Since an obvious way to accelerate 
convergence to a steady state Is to Increase the size of the time step, an 
Implicit scheme Is expected to have a faster convergence. Although the alter- 
nating direction Implicit (ADI) scheme has been valuable In two-dimensional 
problems. Its Inherent limitations In three dimensions suggest an alternative 
approach. A LU Implicit scheme (ref. 13) was demonstrated to be efficient and 
robust for external flows In a broad flow regime and can be readily extended 
from two to three dimensions. This scheme, briefly discussed here, will be 
used for the numerical simulation of turbine flows. 

The conservation law form of the Euler equations In Cartesian coordinates 
for two-dimensional flow Is 


aw 

3F 

3G 

it + 

ix + 

ay 


( 4 ) 


where W Is the vector of dependent variables and F and G are convective 
flux vectors: 

W = (p,pu,pv,pE) T 

2 T 

F = [pu.pu + p,pVU,u(pE + p) 3 (5) 

2 T 

G = [pv , puv , pv + p,v(pE + p) ] 

where p, u, v, E, and p are density, velocity components, total energy, and 
pressure. The pressure Is obtained from the equation of state 

p = p{y - i)[e - 2 + v2 >] < 6 ) 

where y Is the ratio of heat capacity at constant pressure to heat capacity 
at constant volume. These equations are to be solved for a steady-state 
3W/dt = 0, where t denotes time. 


Let the Jacobian matrices be 


a£ aG 

A = aw ’ 8 = aw 

and let the correction be 


6W 

where n denotes the time level. 



The linearized Implicit scheme for a 
equations such as the Euler equations can 


system of nonlinear hyperbolic 
be formulated as 
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( 7 ) 


^1 + B At (D x A + D y B)j AW + At R = 0 
where I Is the Identity matrix and R Is the residual 

R = D x F(W n ) + D y G(W n ) 

Here D x and Dy are central difference operators that approximate 3/3x 
and d/dy. 

If B = 1/2 the scheme remains second-order accurate in time; for other 
values of B the time accuracy drops to first order. The unfactored Implicit 
scheme (eq. (7)) produces a large block-banded matrix, which Is very costly to 
Invert and requires huge storage. An unconditionally stable Implicit scheme 
that has error terms at most of order (At)2 in any number of space dimen- 
sions can be derived by LU factorization 

+ B At (D~A + + D y B + )j [l + & At (D*A~ + D y B") j AW + At R = 0 (8) 

where D* and Dy are backward difference operators and D$ and DJ 
are forward difference operators. The reason for splitting Is to ensure the 
diagonal dominance of lower and upper factors as well as to make use of the 
built-in Implicit dissipation. 

Here A+, A-, B+, and B- are constructed so that the eigenvalues of "+" 
matrices are nonnegative and those of matrices are nonpositive. 

A + = 2 (A + r A I), A = 2 (A - r A I) 

B + = 2 (B + r B I), B' = \ (B - r B I) 

where 

r A i max( |\ A I), r fi i max( |Xg| ) 

Here \a and \b represent eigenvalues of Jacobian matrices. Equation (8) 
can be Inverted In two steps. The LU Implicit scheme needs the Inversion of 
sparse triangular matrices, which can be done efficiently without using large 
storage. This scheme has only two factors In three dimensions. 


CONCLUDING REMARKS 

Two- and three-dimensional composite grids were generated In an attempt to 
Improve grid quality In terms of smoothness, resolution, and orthogonality for 
an annular turbine cascade. This composite grid, which has a C-type (or 
0-type) grid In the Immediate vicinity of the turbine blade, provides a good 
boundary-layer resolution around the leading-edge region for viscous calcula- 
tion, has orthogonality at the blade surface and slope continuity at the C-H 
(or 0-H) Interface, and controls mesh distribution In the upstream region with- 
out using excessive grid points. This composite grid eliminates the undesir- 
able qualities of a single grid when generated for a typical turbine geometry. 

A finite-volume lower-upper (LU) Implicit scheme can be used In solving 
for the turbine flows on the composite grid. This grid has a special grid node 
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that Is connected to more than four neighboring nodes In two dimensions and to 
more than six nodes In three dimensions. But the finite-volume approach poses 
no problem at the special point because each Interior cell has only four 
neighboring cells In two dimensions and only six cells In three dimensions. 

The LU Implicit scheme was proved to be efficient In a broad flow regime and Is 
expected to yield accurate solutions on the Improved composite grid. 
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Figure L - Radial inflow turbine rotor. 




Figure 2. - Surface definition of radial-inflow turbine. 
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Figure 3. - Three-dimensional meshes in radial turbine passage. 







<b) Secondary velocity vectors on cross-channel surfaces. 

Figure 5. - Evolution of passage vortex (from three-dimensional inviscid computation). 







(a) Ring showing cutout for laser window. 
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<b) Geometry at mean section. (V/V cr ■ critical velocity ratio. ) 
Figure 6. - Core turbine stator vane. 
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<b) Two dimensions. 
Figure 7. - Concluded. 
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(a I Two dimensions. 

Figure 8. - Composite grid for turbine stator vanes of annular cascade. 
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(b) Cell centers to be used in finite-volume scheme. 
Figure 10, - Special point of composite grid. 
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